Division of the two-qubit Hilbert space according to the entanglement sudden death 

under composite noise environment 
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We show theoretically that according to the disentanglement behavior under composite noise 
environment, the Hilbert space of a two-qubit system can be divided into two separate parts: a 3- 
dimensional subspace in which all states disentangle asymptotically, and the rest in which all states 
disentangle abruptly. The violation of additivity for entanglement decay rates under weak noises 
[see, PRL 97, 140403 (2006)] therefore can be explained in terms of such division of the Hilbert 
space. 
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I. INTRODUCTION 



Quantum coherence, as a consequence of the super- 
position principle, lies in the heart of quantum mechan- 
ics and is regarded as a unique sign for the quantum 
regime. A fundamental question raised by quantum co- 
herence is the the emergence of macroscopic classical- 
ity from the microscopic quantum world. It is currently 
well-recognized that the answer lies in the so-called de- 
coherence process, a dynamical evolution of the system 
in which quantum coherence is gradually lost due to 
the ubiquitous system-environment interaction De- 
coherence has been extensively studied for years in the 
realm of quantum optics as a decay of local co- 
herence. However, it was not until recently that the 
evolution of nonlocal quantum coherence, i.e., the en- 
tanglement dynamics, has become a field of interest 
that attracts growin g th eoretical and experimental stud- 
ies SiiiaaiMtiltiltilQEKtiS- More- 
over, since entanglement itself is a key resource for quan- 
tum computation and information processing [TR [l9| . 
such studies not only put new insights to the funda- 
mentals of quantum mechanics, but also scrutinize more 
stringently on whether and how one can build an appli- 
cable quantum computer under realistic circumstances. 

A remarkable finding in the study of entanglement dy- 
namics is the entanglement sudden death (ESD), which 
was previously shown in [3j and latterly coined by Yu and 
Eberly [H0]. The model studied in Q consists of two 
initially entangled but non-interacting qubits each cou- 
pled to its own local environment. In all circumstances 
the entanglement between the two qubits disappears sub- 
sequently due to the spontaneous decay of each qubit. 
Very interestingly, the entanglement dynamics falls into 
two distinct categories: abrupt and asymptotic disentan- 
glement, depending on the initial states. The abrupt 
disentanglement is characterized by a complete disap- 
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pearance of entanglement within finite time and is thus 
termed ESD. The significance of ESD lies in two- folds: 
(i) It unveils a fundamental difference between local and 
nonlocal quantum coherences via their evolution dynam- 
ics; (ii) It puts an upper bound on the applicability of en- 
tangled pairs in practical quantum communication even 
with the best protocol for entanglement distillation. 

Besides the apparent dependence on initial states, en- 
tanglement dynamics is also intrinsically impacted by the 
specific form of the system-environment interaction. Yu 
and Eberly also demonstrated examples of ESD under 
composite noise environment [5| and classical noise envi- 
ronment [f|. Structured reservoirs with memory effects 
add a new ingredient, entanglement revival, to the entan- 
glement dynamics, due to the non-Markovian nature of 
the system- reservoir interation [7j. Ikram et al studied 
the effects of squeezed and thermal reservior Q. Al- 
Qasimi et al @ showed that all X-states undergo ESD 
when the reservoir is at finite temperature. Extensions 
of Yu and Eberly's model to commonly shared environ- 
ment (Dicke- regime) (lpj and high dimensional bipartite 
systems (e.g., the 3 <£> 3 system [ll[ and the continu- 
ous variable system [l2j]) have also been studied recently. 
Moreover, the question on where the lost entanglement 
goes has been addressed recently [H, [l4| , revealing a rich 
and counterintuitive relation between the ESD and the 
entanglement sudden birth in the reservoirs. 

However, in most studies 0, S 1, 0, S, 1 the two-qubit 
state under investigation is artificially confined to certain 
simple classes (e.g., Bell- like states, Werner-like states, 
or the X-states) for the convenience of theoretical han- 
dling. This led to a lack of global knowledge on how the 
full 4-dimensional two-qubit Hilbert space is structured 
according to the two distinct behaviors of entanglement 
decay. From the perspective of quantum communication, 
it is very important to know which state is more robust 
than others in fighting against entanglement lost. Thus it 
is desirable to discriminate states with asymptotic disen- 
tanglement behavior from others. Very recently, Huang 
and Zhu [l5[ have specified the necessary and sufficient 
condition for ESD under amplitude damping and phase 
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damping, respectively, using a principal-minors technique 
based on the positive partial transpose (PPT) criterion. 
The purpose of this paper is to show that, even under 
the combined action of these two noises, there exists a 
3-dimensional ESD-free subspace, in which all state dis- 
entangles asymptotically. Furthermore, we have proved 
that all states outside this subspace disentangle abruptly 
under the composite noise environment. Thus the full 
Hilbert space is completely divided into two separate 
parts determined solely upon whether they undergo ESD 
or not. Such a division of the Hilbert space also provides 
a simple explanation to the violation of additivity for 
entanglement decay rates under weak noises @ via the 
restriction of the ESD-free subspace on the system dy- 
namics. 

This paper is organized as follows: In Section |TT] we 
prove that the combined action of amplitude damping 
and phase damping can be factored as subsequent ac- 
tions of each individually. Therefore the set of ESD-free 
states under composite noise environment is contained 
in the intersection of the two sets of ESD-free states un- 
der each of the constituent noises. In Section imi we show 
that under pure phase damping, the ESD-free states con- 
stitute four 3-dimensional subspaces. In Section IIVI we 
prove that only one of the four subspaces is ESD-free un- 
der the composite noise environment. Finally, we make 
a few remarks and give a short conclusion in Section fVl 



II. THE FACTORIZATION OF AMPLITUDE 
DAMPING AND PHASE DAMPING 

The system we study includes two non-interacting 
qubits (A and B) and two independent local environ- 
ments, which will be modeled as either amplitude damp- 
ing, or phase damping, or the combined action of both. 

The master-equation of the system under amplitude 
damping and phase damping can be written as 



operators as 



p(t)= ]T K iP (0)Kj, 



(3) 



i=l,2,3,4 



with = 1- The two-qubit Kraus-operator 

KiS can be expressed as tensor product of the one- 
qubit Kraus-operators for each party. We specifically 
denote Kf M as the two-qubit Kraus-operator under am- 
plitude damping defined by Kf m = ® Aff, K^ m = 



Mf ® Mf , Kg 
M^ B = 
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where ^f :B = exp(— ^Tf' B t) and u^ B = y 1 — {jf' B ) 2 - 



While the Kraus-operator K[ for phase damping is de- 
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where t^' S 



exp(— T 2 ' B t) and u) 2 ' E 
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When both noises participate, the master-equation reads 
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Similarly, the time evolution of Eq. ( [6]) can be expressed 
by the Kraus-operators as 

p(t)= E (c£®cf) P (o)(c£®cpy , (7) 



where 
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respectively. The two-qubit density operator p is repre- 
sented in the computational basis |1) = | ||), |2) = | fj,), 
|3) = | It), |4) = | U). The parameters V( ' B and T^ B 
(for party A, B) are the population relaxation rate and 
dephasing rate for amplitude damping and phase damp- 
ing, respectively. The operators cr+' B , a^' B and <r^' B 
denote the raising, lowering and Pauli operators for each 
qubit. The explicit time evolution of Eq.([T]) and can 
be solved and expressed in a unified form by the Kraus- 
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Using the above Kraus-operators it is straightforward 
to verify the following equations 

£ (c£®cf)p(o)(c£®cfy 
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which can be further packaged into a compact form as 



(10) 



where $y m , $f , and $f denote the quantum channels 



governing the evolution of the system under amplitude 
damping, phase damping, and both of them, respectively. 
The essence of Eq . ( fTUj) is that the effect of the compos- 
ite noise environment can be factored as subsequent ac- 
tions of each constituent noise environment. Moreover, 
as a consequence of the Markovian master-equations, it 
is also straightforward to verify a general time-domain 
factorization property bellow 

$ t = $ T /$ r , (11) 

where t = t + t', t,t' > 0. Eq.( fTTj) holds for every 
kind of the quantum channels $f , $f m and $f . Since 
the entanglement is non-increasing under local operations 
(note that all the quoted quantum channels are local op- 
erations) [HI, [lj|, a direct conclusion drawn from the 
factorization properties Eq.( [TO]) and ( [TT|) is that if a 
state has completely disentangled under either %f- m or 
$f l , it must remain disentangled under the action of $f . 
In other words, the set of the ESD-free states under the 
composite noise environment is a subset of the intersec- 
tion of the sets of ESD-free states under each one of the 
noises. Therefore, to be ESD-free under the action of 
each noise is a necessary condition to be ESD-free under 
the action of both noises. 

Finally, we would like to clarify that amplitude damp- 
ing itself does not merely produce damping in the ex- 
cited state of the qubit, but also causes dephasing be- 
tween the two states. This concomitant dephasing has 
been automatically included in the the master-equation 
([T|) and the quantum channel %f m . In other words, it is 
impossible to factor $f m further into a fictitious "pure 
amplitude damping" channel and a fictitious "pure de- 
phasing" channel. Therefore, when we talk about the 
composite noise environment we mean that an extra de- 
phasing noise, besides the dephasing caused by amplitude 
damping, is present in the system-reservoir interaction. 
Such composite noise environment is a reasonable model 
for realistic circumstances where the noisy environment 
can be eventually treated as a mixture of the two noises. 



III. THE ESD-FREE SUBSPACE UNDER 
PHASE DAMPING 

We have shown that to find the ESD-free states un- 
der the composite noise environment, it is sufficient to 
explore within the set of the ESD-free states under ei- 
ther one of the two noises. Our strategy is to investigate 
the set of ESD-free states under phase damping first, 
since there is a simple partition between the ESD and 
the ESD-free states for pure phase noise. The question 
for the composite noise environment will be addressed in 
the next section. 



The necessary and sufficient condition for ESD un- 
der phase damping has been specified by Huang and 
Zhu [l5[ very recently using the PPT criterion. The 
method they use is technically tedious and lack of phys- 
ical transparency. For the completeness of narration, we 
re-derive it u sing a more intelligible method based on the 



concurrence 



20l | . The concurrence C{p) is an entangle- 



ment monotone ranges from to 1 and can be computed 
directly from the two-qubit density matrix p as 



C{p) = max(0, A), 



(12) 



where A = \f\\ — \f\~2 — V^3 — V^4> with the Aj's the 
eigenvalues of the matrix 



R{p) = p (a y (g>(Ty)p* (dy ® a y ) 



(13) 



in descending order. 

We first consider the concurrence under phase damping 
in the infinite-time limit, i.e., under the action of 
Suppose the initial state is 



p(0) = 



Pll Pl2 Pl3 Pli 

P21 Pl2 P23 P24 

P31 P32 P33 P34 

P41 P42 P43 P44 



(14) 



After infinite time of phase damping it is found 
p(oo) = $™p(0) = 



pn 

P 22 

p 33 

P44 



The concurrence of the above final state is simply 
C[p(oo)] = max [0, A(oo)] , 



where 



A(oo) 



-2y/pUpU, P22P33 > P11P445 



-2^/522/033, P22P33 < P11P44- 

Suppose P11P22P33P44 7^ 0, it is obvious that 
A(oo) < 0. 



(15) 



(16) 



(17) 



(18) 



Since A(t) is an algebraic function of the eigenvalues of 
R[p(t)], it must be a continuous function of t. The in- 
equality (fl~8|) thus means that if A(t) starts from a posi- 
tive value, it must have crossed over zero before t = 00, 
i.e., the two-qubit system must have disentangled at a 
finite time. That is to say, the necessary condition for a 
state to be ESD-free under phase damping is to have at 
least one vanishing diagonal element in its density ma- 
trix. Thus any ESD-free state must belong to one of the 
four subspaces below 



Pi = 





P22 P23 P24 

P32 P33 P34 

P42 P43 P44 
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Pin 



Piv 



' Pll p 13 pu 



P31 p 33 p 34 

. PAl p 43 p 44 

Pll Pl2 pu 

P21 P22 ,024 



Pll P42 p 44 

Pll Pl2 Pl3 ' 

P21 P22 P23 

P31 P32 P33 





(19) 



To prove that all states in pi^iv are indeed ESD- 
free under phase damping, we further investigate the 
concurrence at finite time t, e.g., for an initial state 
pi(0) located in the subspace pj. We find the eigen- 
values of R[pi(t)\ = i?[$fV/(0)] are Ai, 2 = a ± (3, 
A3, 4 = 0, where a = {72I2 ) 2 P23P32 + P22P33 and = 

2 l2l2 \J P23P32P22P33- ThuS 



(20) 



A(i) = VAi - VA 2 - VA 3 - V A 4 
= 2j3/Ma + P+y/a-p). 



It is clear that A(i) > holds Vt. Only in two cases 
A(t) — 0: (i) t — 00, which means asymptotic disen- 
tanglement; (ii) P23 — P32 — 0, which means A(t) = 
holds Vi, i.e., the state pi{0) is seperable from begin- 
ning. Thus all states in pi are indeed ESD-free under 
phase damping. Straightforward calculations show that 
the same conclusion holds for the other three subspaces 
Pi 1, Pi 1 1 and piy. Thus, a generic state p(0) is ESD-free 
under phase damping iff 



P11P22P33/O44 = 0, 



(21) 



or equivalent to say, iff it is located within the four sub- 
spaces pi^iv shown in Eq.f [H 



IV. THE ESD-FREE SUBSPACE UNDER 
COMPOSITE NOISE ENVIRONMENT 

Now we further explore the entanglement dynamics for 
states in the subspace pi^iv under the combined action 
of phase damping and amplitude damping. We find com- 
pletely different disentanglement behaviors between the 
subspace pi and the other three subspaces pn^iv'- all 
states in pi disentangle asymptotically and all states in 
Pi 1, Pi 1 1 and piv disentangle abruptly. 



A. All states in pi disentangle asymptotically 

It is straightforward to show that the subspace pi is 
full ESD-free whatever the noisy environment is. For 



amplitude damping, the evolution of an arbitrary ini- 
tial state pi(0) is expressed as p' 7 (i) = $^ m yO/(0) and 
the eigenvalues of i?[p' 7 (i)] are A^ 2 = a' ± and 



a; 



3,4 



0, where a' 



(7i 7i ) (P23P32 + P22P33) and 
0' = 2(7^7f ) 2 ^/P23P32P22P33- Obviously, (3' > and 
the equality holds only for t = 00 or P23 = P32 = 0. 
Following the same reasoning as for Eq.f [20]) . it is clear 
that the subspace pi is also ESD-free under amplitude 
damping. Similarly, when both noises are present, the 
evoluation of pi(0) is simply p"{t) = $fpj(0) and the 
eigenvalues of R[p"(t)] are A" 2 = a" ± [3" and A 3 ' 4 = 0, 
where a" = (7j Sf ) 2 [(l2li) 2 P23P32 + P22P33] and f3" = 
2hilF) 2 l2li\/P23P32P22P33- Again, (3" > and the 
equality holds iff t = 00 or p 2 3 — P32 = 0. Thus pi is also 
ESD-free under composite noise environment. 



B. All states in pn, pin and piv disentangle 
abruptly 

The situations in pn^iv are more complicated because 
the explicit expression of concurrence at finite time t is 
too lengthy to give a simple conclusion (due to the fact 
that the forms of pn^iv do not preserve under ampli- 
tude damping). Such complexity is also reflected in [15l |. 
where the derived necessary and sufficient condition for 
ESD under amplitude damping is mathematically pack- 
aged into the question of the positive definiteness for an 
artificial 4x4 matrix {see Eq.(ll) of [H|]}, which is still 
far from transparent. 

However, using the factorization properties we can in- 
deed prove that no state in pn^iv is ESD-free under 
composite noise environment. The trick is to divide the 
quantum channel $p into two subsequent quantum chan- 
nels $^ and as in Eq.( [TT|) . The evolution of a generic 
intial state p(0) under 



thus abides 



p(t) = $^p(0) = $^p(0) = $>(r), (22) 

where p(r) — $^p(0). Therefore the final state p(t) can 
be obtained by acting $ T / over the intermediate state 
p{r). Now one asks what is the necessary and sufficient 
condition for asymptotic disentanglement for p(0) under 
the action of $p? The answer is quite simple: Neither 
has p(r) completely disentangled within finite time t , nor 
shall p[r) disentangle completely under the action of$^,, 
Vt, t' > 0. The first part of the statement is no more 
than a repeat of the original question. However, the sec- 
ond part of the statement requires that over the full time 
evolution any intermediate state p(r) should belong to 
the set of the ESD-free states under the action of $^,. 
Since the set of the ESD-free states under S^", is contained 
in the subspaces pi^iv, it is necessary for any interme- 
diate state p(r) to be located within pi^iv to guarantee 
asymptotic disentanglement. Now it is straightforward 
to show that as long as p\\ 7^ 0, such requirement can- 
not be satisfied whenever amplitude damping is present 
for each party: the diagonal elements of p(r) are solely 
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determined by amplitude damping as 

Pii(r) = (7iV) 2 Pii, 

p 22 (r) = (7i A ) 2 [P22 + H B ) 2 PiiL 

P3 3 (r) = (7f) 2 [P33 + (^) 2 pn], 

PAa{t) = PAA + (wf ) 2 P33 + (w 1 1 ) 2 P22 + (wfwf ) 2 pn, 

(23) 

where 7^"' S , are functions of r instead of t. It is 

obvious that if pn ^ and I\ ' ^ 0, none of the above 
diagonal elements vanishes Vt > 0. Therefore, we arrive 
at the conclusion that no state in pn^iv is ESD-free due 
to the existence of nonzero element p\\ (Note that there 
exist states with pn — in pn^iv, however, this trivial 
ambiguity can be excluded by absorbing all these states 
into pi). 

V. DISCUSSION AND CONCLUSION 

An interesting topic connected with our work is the 
violation of additivity for entanglement decay rates un- 
der weak noises discovered by Yu and Eberly They 
have demonstrated that a state disentangles asymptoti- 
cally under either amplitude damping or phase damping 
may, however, disentangles abruptly under the combined 
action of both. This fact unequivocally manifests the vi- 
olation of additivity for the decay rates of nonlocal quan- 
tum coherence, in sharp contrast to the uphold of that 
for local quantum coherences. In view of the ESD-free 
subspaces discussed, such phenomenon deserves a sim- 
ple explanation. The specific example used in [f| is no 
more than a special state located in the subspace piv- 
Although it is ESD-free under either one of the noises, it 
is not ESD-free under both because the evolution of the 
state (due to amplitude damping) is not confined in the 
ESD-free subspace required by the composite noise envi- 
ronment. According to the previous analysis we conclude 
that such an abrupt violation of additivity can only occur 



for states located in the subspaces pn^iv ■ On the other 
hand, due to the existence of the completely ESD-free 
subspace pi, it is still possible to explore some quasi- 
additivity of entanglement decay rates for states within 
this restricted space. 

Another important work to note is the finite tem- 
perature effect recently considered by Al-Qasimi and 
James @ for amplitude damping (i.e., spontaneous exci- 
tation from | |) to | t) is also allowed due to thermal exci- 
tations in the reservoir). As demonstrated in an X-state 
example, they argue that all states disentangle asymp- 
totically in zero-temperature bath must undergo sudden 
death at finite temperature. This is a very strong state- 
ment that may smear off the necessity to discriminate 
between the two kinds of disentanglement behaviors un- 
der realistic circumstances. Therefore, the significance 
of our work is restricted in the zero-temperature regime, 
where spontaneous excitation is forbidden. 

In conclusion, we have shown that the set of ESD-free 
states for a two-qubit system under the combined action 
of amplitude damping and phase damping is fully repre- 
sented by the subspace pj, which is spanned on the bases 
I TO; I IT) an d I 10- Therefore, any contamination from 
the double-excitation state | |T) i n the initial density ma- 
trix is extremely hazardous for entanglement preserva- 
tion. Entanglement resources realized via | ft) ( e -g-; the 
Bell-state <fi±) thus should be treated with great caution 
comparing to those without it (e.g., the Bell-states tjj±) 
whenever the lifetime of entanglement has to be taken 
into consideration. 
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